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Abstract 



We derive the semiclassical series for the partition function in Quantum 
Statistical Mechanics (QSM) from its path integral representation. Each term 
of the series is obtained explicitly from the (real) minima of the classical 
action. The method yields a simple derivation of the exact result for the 
harmonic oscillator, and an accurate estimate of ground-state energy and 
specific heat for a single-well quartic anharmonic oscillator. As QSM can be 
regarded as finite temperature field theory at a point, we make use of the field- 
theoretic language of Feynman diagrams to illustrate the non-perturbative 
character of the series: it contains all powers of h~ and graphs with any number 
of loops; the usual perturbative series corresponds to a subset of the diagrams 
of the semiclassical series. We comment on the application of our results 
to other potentials, to correlation functions and to field theories in higher 
dimensions. 
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I. INTRODUCTION 



Semiclassical series have a long history in Quantum Mechanics which goes back to the 
early days of the Schrodinger equation. Their first terms are dictated by classical trajectories, 
and serve as the initial step of iterative procedures which yield all other terms [ij. 

Path integral representations for correlation functions have been developed more recently 
@-||. Being sums over trajectories, they provide a natural derivation of semiclassical results 
through the stationary phase method. Indeed, the saddle points of the path integrals cor- 
respond to classical trajectories which dictate the first terms of the series However, 
despite the many applications in Quantum Mechanics and Field Theory, most discussions 
which used path integrals never went beyond the first term of a semiclassical series. Notable 
exceptions were the works of DeWitt-Morette |]ID| and Mizrahi |TT| in Quantum Mechanics. 

Semiclassical methods for finite temperature field theories [|12]-[14j also remained re- 
stricted to derivations of the first term of a semiclassical series [fL5|J , even when the problem 
was reduced to Quantum Statistical Mechanics fll||r7| , viewed as field theory at a point (zero 
spatial dimension). Some references resorted to extensions to the complex plane |T8|-pO| to 
include complex paths required to describe Fourier transformed quantities but, again, those 
treatments were not concerned with obtaining the whole series. 

In this talk we will present a systematic path integral procedure to generate semiclassical 
series in Quantum Statistical Mechanics. It leads to the construction of each term of the 
series from the solution(s) of the classical equations of motion. We will focus our attention 
on the partition function, and use the method of steepest descent which, in this case, only 
requires real solutions as saddle-points W% - The restriction of our analysis to quantum- 
mechanical systems (i.e., field theories at a point and at finite temperature) will allow us to 
construct the semiclassical propagator needed to generate the terms of the series. 

In our contribution to last year's Workshop ||22|| , we had already outlined the procedure 
mentioned in the previous paragraph. However, the approach we used then was quantum- 
mechanical, as opposed to the field-theoretic language we shall adopt in the present article. 
Reference describes both approaches and gives a detailed account of the results which 
will be quoted here. 

This article is organized as follows. Section |TJ presents the derivation of the semiclas- 
sical series for a generic one-dimensional potential of the single-well type in field-theoretic 



language, which allows for a simple connection with the works of references [10,11|; the 



presentation is a natural extension of textbook material [23], and profits from the clear 



account of reference P0| . Section [TTT] applies our results to the harmonic oscillator and to 
the single-well quartic anharmonic oscillator; for the latter, we compute the ground-state 
energy and the specific heat. Section [TV] presents our conclusions, comments on extensions 
to double-well type potentials, to QSM in higher dimensions and to field theories. 



II. QUANTUM STATISTICAL MECHANICS 

The partition function for a one-dimensional quantum-mechanical system consisting of a 
particle of mass m in the presence of a potential V(x) in equilibrium with a thermal reservoir 
at temperature (3~ l can be written as a path integral: 
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30 rx(/3%)=XQ 

dx / [Vx(t)\ e~ s / h , 

oo Jx(0)=xq 



(1) 



S[x) 



dr 



o 



HI) +v{x) 



(2) 



For convenience we define the dimensionless quantities q 

U(q) = V(x]siq) /muj 2 N x 2 N and g = h/mu N x%, where u^ 1 and x^ are natural time and length 
scales of the problem, respectively, 
function as 



x/x N , 6 = UJ N T, = /3Hun, 
• N are na 

In terms of these quantities we rewrite the partition 



Z(Q) 



dq / [Vq{6)) e"^, 

Jq(0)=q 



(3) 



e 



dB 



^Q 2 + U(q) 



(4) 



where the dot denotes differentiation in 9. 

We generate a semiclassical series for Z(Q) by: (i) finding the minima q c (6) of the 
Euclidean action /, i.e., the stable classical paths that solve the Euler-Lagrange equation 
of motion, subject to the boundary conditions; (ii) expanding the Euclidean action around 
these classical paths; (iii) deriving a quadratic semiclassical propagator by neglecting terms 
higher than second order in the expansion; (iv) using that propagator to compute higher 
(than quadratic) order contributions perturbatively. 

For the sake of simplicity, we shall restrict our analysis to potentials of the single-well 
type, twice differentiable, and such that U'(q) = only at the minimum of U, which we shall 
assume to be at the origin (see Fig. [I]). This guarantees that, given go and 0, there will be 
a unique classical path satisfying the boundary conditions. Multiple-well potentials force us 
to consider more than one classical path for certain choices of go and 0. This phenomenon 
has been analyzed, for a double-well type potential, using the language of catastrophes and 
bifurcations 



Semiclassical series for the double-well quartic oscillator will be presented 
The Euler-Lagrange equation (U' = dU/dq) 



elsewhere EI 



q-U'(q)=0, 



(5) 



subject to the boundary conditions q(0) = q(Q) = qo, describes the motion of a particle in 
the potential minus U. Its first integral is 



q 2 = U(q)-U(q t ), 



(6) 



where q% denotes the single turning point (since we have an inverted single well) of the 
motion, defined implicitly by 



= 2 



dq 



in 



v(q,qt)' 



(7) 
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where v(q,q') = sign(g' — q)^2[U(q) — U(q')]; equation (|7|) is a consequence of integrating 
(§). Thus, for a single well, given go and 0, the classical path will go from go, at 9 = 0, to 
It = qt(qo, @), at = 0/2, and return to go at = 0. (Note that sign(g 4 ) = sign(g ).) 
The action for this classical path has a simple expression in terms of its turning point: 

rqt 

I[q c ] = eU(q t )+2 dqv(q,q t ), (8) 
Jqo 

where we have used (|6|). The first term in (|8[) corresponds to the high-temperature limit of 
Z(Q), where classical paths collapse to a point (q t — > g ). The last term will be negligible for 



potentials that vary little over a thermal wavelength A = hy (3/m. However, by decreasing 
the temperature it will become important and bring in quantum effects. 

We now expand the action around the classical path. Letting q{9) = q c {9) + r](9), with 
77(0) = r?(0) = 0, we obtain 



where 



I[q]=I[q c )+I 2 [r)} + 6I[r)}, 



hH = \[de {rf(0) + U"[q c (9)]v 2 (0)} 



(9) 



(10) 



5I[rj] = f d95U(9 lV ) = £ - / d9U^[q c (9)}r ] n (9). 
Jo n=3 n. Jo 

Inserting (^) into (|3|) and expanding e~ 5I ^ g in a power series yields 
The summation in ( |12|) can be written more explicitly as 



£ 

m=0 



si[ V y 



i+£ 



n 



00 , e 

£ — L d9 3 U^[q c (9 j )] V ^(9 J 



n,=3 n i ! J0 



9 ) ^1 9"*m\ "1 

As a consequence, one is led to compute integrals of the following type: 

/■i?(e)=o 

( V (9 1 )--- V (9 k )) ee / [P7 7 (e)]e- /2 [°^]/^(e 1 ).--7 7 (^). 

Jr;(0)=0 



(11) 



(12) 



(13) 



(14) 



Such integrals emerge naturally as functional derivatives of the following generating func- 
tional: 



Z[J] 



r,(6)=0_ 
r,(0)=0 



[V V (9)} e 



•i|j a [0,e W ]-/ o e d9J(tf)i7(tf)} 



(15) 



Indeed, 
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(r ] (9 1 )...r ] (9 k ))=g k 



In order to compute 2[J], we define 



5 k Z[J] 



5J(9 1 )---5J{6 k ) 



j=o 



e 



v (9)=fj(9) + / dO'g{6,0')J(6') 
Jo 



where r/(0) = 77(0) = 0, and G(9,9 f ) satisfies 

r ^ 



a^ 2 



tf"[?c(0)] G(9, 9') = 5(9 - 9'), £(0, 9') = £(0, 9') = 0. 



Inserting ( |17|) in (|15[) , and noting that [X>?y(0)] = [X>?y(0)], we obtain 



/?(e)=o 



Jij(0)=0 



If we define 



G c (9 1 ,r] 1 ;9 2 ,r]2) 



r)(02)=ri2 , 
r?(0l)=m 



-h[6i,6 2 ;r,]/9 



(16) 



(17) 



(18) 



(19) 



(20) 



h[Bx^n] = -J d9 {n 2 + u"[q c (9)]v 2 }, 



we finally arrive at 

Z[J] = G c (0,0;9,0) exp 

Using this result, we can now calculate (|T6"D. The result is simply 

(V(0i) ■ ■ ■ V(9 k )) = g k/2 G c (0, 0; 6, 0) £ G(9 h A 2 ) ■ ■ ■ G(9 lk _ 17 9 lk ), 



— [ d9 [ d9'J(9)g(9,9')J(9') 
2g Jo Jo 



(21) 



(22) 



(23) 



if k is even, and zero otherwise. J2p denotes sum over all possible pairings of the 9i .. 
Inserting this into (|P2|) and (|H|) yields the semi classical series for Z(Q). 

We still have to solve Eq. (JT^). This can be easily done if one notes that, for 9 7^ 9', it is a 
homogeneous second-order differential equation. Therefore, G(9, 9') can be constructed from 
a linear combination of two linearly independent solutions rj a (9) and r]b(9) of the equation 



Indeed 



ii-U"[q c (9)]v = 0- 



nia a>\ - S a-Va(9) + b^ Vb (9), 9 < 9' 
y[ ^ U) \a +Va (9) + b +Vb (9), 9>9'. 



(24) 



(25) 



Continuity imposes 
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Q{Q> + e,9') = G{& -e,9'), 



(26) 



whereas 



leads to 



d_ 
36 



d_ 

06 



0(0,0'] 



(27) 



with e — > + . (Pq), ( p7|) and the boundary conditions completely determine the coefficients 
in (|25|). The final result is 



Q(0,0') 



Q(o,e < )n(e > ,Q) 



o(o, e) 

where 0<(0>) = min(max){6 1 , #'}, and f2(#i,#2) is the function 

n{6,e')=r] a (6) Vt^-Vam^d). 



(28) 



(29) 



In the Appendix we show that G c (6i, rji, 62, ^2) can also be obtained from the two linearly 
independent solutions of Eq. (0), f] a {6) and rjb(9)] furthermore, we show how to construct 
those two functions from the solution q c (9) of the classical equation of motion. This com- 
pletes the steps needed to write down any term of the series: all that is required is q c {0)\ 



III. QUANTUM OSCILLATORS 

In this section we will apply our construction to the harmonic oscillator and to the single- 
well quartic oscillator. The harmonic case is designed to illustrate the compactness of our 
general formulae, which immediately yield the exact answer — there is no need to compute 
functional determinants from eigenvalue problems! The anharmonic case is designed to 
illustrate their power — the first term of the semiclassical series for the partition function 
allows us to extract a very good estimate of the ground-state energy and of the specific heat. 



A. The Harmonic Oscillator 



In this subsection, we study the potential 

1 



V(x) = ^mu 2 x 2 . (30) 



Choosing uo N = uo and xn = yh/mu, and introducing the dimensionless quantities of section 
U we have g = 1 and 



Integrating (|) leads to 



U(q) = \q 2 . (31) 



q c (6) = q t cosh(0-|), (32) 
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The relation between qo and qt is obtained by taking 9 = in (|32|): 

<?o = g c (Q) = % cosh(6/2). 
The action for the classical path is 

I[q c ] = g o 2 tanh(0/2). 
Following the Appendix, the functions rj a and rjf, are given by 

e, 



and 



Va(6) = q c (9) = Qt sinh(6> 



Vb{0) = qje) QW) = 'Qt 1 cosh(0 - 



e, 



(33) 



(34) 



(35) 



(36) 



It follows that = sinh(#j — Qj) and = cosh(9j — Oi), leading to 

1 



G c (9x,r]x; 6*2,772) 



27Tsinh(^ 2 - Ox 



x exp 



1 



2sinh(fl 2 - 9x) 



cosh(6» 2 - 9 X ) {rjl + Vi) ~ 2 V1V2 



(37) 



Using (|A10|) we obtain A = 2n sinh6. Since the problem is quadratic, its exact solution is 
then given by 



/oo 
dq e ~nic]/g^-i/2 
-00 

Inserting ( 34|) and the value of A yields the well-known result 



(38) 



Z(B) 



-ql tanh(e/2) 



dq 



V2n sinh 2sinh(6/2)' 



(39) 



B. The Single-well Quartic Oscillator 



In this subsection, we study the potential 

1 



V(x) = — muj 2 x 2 + - Ax 4 . 



(40) 



Choosing = ui and xjy = \Jmtu 2 /X, and introducing the dimensionless quantities of 
section ||, we have g = \h/m 2 u) 3 and 



(41) 



Integrating (^) leads to [f25| , f26 



Qc(0) = qtnc(u e , k), 



where nc(u, k) = l/cn(w, k) is one of the Jacobian Elliptic functions [25-^7|, and 

0' 



u e = Jl + qi[9 



k 



2 + qt 



2j' \ni+q?y 

For future use, we note that (|43|) can be rewritten as 

29 -0 



(42) 



(43) 



2V2P-1' 



VIA 



'2(1 -k 2 ) 
2k 2 - 1 ' 



The relation between g an d % is obtained by taking 9 = in (32): 

?o = ?c(©) = 9tnc« e . 

(We shall often omit the /c-dependence in the Jacobian Elliptic functions. 
The action for the classical path (|^) is 



(44) 



(45) 



I[q c ] = I7(<ft) + V2 [ M dq ^^ + q 2 + 2)(q 

J\qt\ 



' 2 - Qt) 



(46) 



Performing the integral (Ref. f25| , formula 3.155.6) and replacing go by the r.h.s. of (f35|), we 
obtain 



+ sn w e ( 1 + - q 2 t nc 2 w e ) y 1 + ^ Qt (1 + nc 2 w e ) 



(47) 



where E(<p, fc) denotes the Elliptic Integral of the Second Kind and (pg = axccos[q c (9) / q ] = 
arccos(cnu ). 

For the construction of the quadratic semiclassical propagators G c (9, rj; 9', rf) and Q(9, 9') 
we shall need 



Va{9) = q c {9) = qtJl + q 2 snu e dnu s nc ug 



(48) 



and 



Q(9) = C 2 (l + % 2 )- 3/2 



l-^)t*+(^-2)E( Wl *) 



cn ue dn ,, 2 ^ cn uq sn u e 



snwe 



+ (k 2 ~ I) 



dnu g 



(49) 



We may then obtain rj b (9) = q c {9) Q{9) and, thus, f2 12 and W\ 2 from ( |X6| ) and ( |^7|) . Finally, 
use of (|A9|) and (E3) will yield the desired propagators. 
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For the series expansion of the partition function, we shall need 



5U(6, V ) = q c (6) V 3 + ^ V \ (50) 

obtained from ( JTT|) . Therefore, we have to consider not only the usual quartic vertex, but 
an additional time(#)-dependent cubic term. This completes the set of ingredients needed 
to write down a semiclassical series for any correlation function. In the next subsection, 
we shall concentrate on the first term of the series for Z(Q), which yields the quadratic 
approximation. 

1. The quadratic approximation for Z(Q) 
From the knowledge of the classical action and of the Van Vleck determinant, we define 

/oo 
dg e- 7 ^ A -i/2 (51) 
-oo 

as the quadratic approximation to Z(Q). To perform the integral over go one must write 
I[q c ] and A solely in terms of go (and 9), but except in rare cases this is not an easy task. 
Usually, it is much simpler to write these quantities in terms of q t [see Eq. (pE5|)], and so it is 
natural to trade go for qt as the integration variable in (|51|). This is much simplified by the 
fact that the Jacobian of the map g — > qt is simply related to the van Vleck determinant. 
In fact, Eqs. (0) and QMT]) imply 

= \v =-T^A (52) 

\dq t J e (dO/dq ) gt 2 \ d 1tJ qo 4:7rgv(q ,q t ) 

Eq. fl5Ip then becomes 

Z 2 (e) = - J- [ q+@ dq t U ' { f Al / 2 e- 1 ^ = [ q+& dq t D(q t , 6) e "^/» (53) 
^9 he v{q ,q t ) J q ~ 

where g| = lim g0 ^ ±oo g t (g , 6). 

The expression above is valid for single-well potentials in general. Now, let us especialize 
to the potential flU"]) . I[q c ] is given by (|47|) , and using (|A10|) and fl49|) one can write D(q t ,Q) 
as 



— «e + — p — ^(Vg, «) 



cnwednwe ,„ , 9 . cn-u« sn-uo 

H + (1 - fc ) 

sn Me dn uq 



1/2 



(54) 



From ( f4"5D it follows that go — > oo when cn(ue, fc) = 0, which occurs when tie = K(/c), where 
K(/c) is the Complete Elliptic Integral of the First Kind. Using (f44|), this condition can be 
written as an equation in k: 



9 



0_ 



K(k). 



(55) 



The graph of f(k) = 2\j2k 2 — 1 K(k) is plotted in Fig. |2|. It increases monotonically from 
zero (at k = 1/V2) to infinity (as k — > 1), and so for each nonnegative value of O Eq. (|55"D 
has a unique solution, which we denote by k®. Eq. (f44|) then gives the corresponding value 
of q£ (q e = -q£, since U(-q) = U(q)). 



2. Limiting cases of the quadratic approximation 



Expression ([53]) may be used to compute Z 2 (Q) numerically for any value of 6. How- 
ever, certain limiting cases may be dealt with analytically. These limits are: the harmonic 
oscillator (g — > 0), high temperatures (6 — > 0), and low temperatures (6 — > oo). 

The limit g — > of (|53|) does yield the partition function of the harmonic oscillator, as 
required, since V(x) = \ muj 1 x 1 when g = 0. In order to arrive at this result, we note that 
in this limit one can perform the integral ([5^) using the steepest descent method. Details 
of the derivation can be found in [|l|]. 

At high temperatures, — > and ( |55| ) is solved for fc — > l/v2, and so g@ — > oo. It 
follows that 



Z 2 (Q) 9 ~° 



1 



2tt#0 



dq e 



-&U(q)/g 



(56) 



or, equivalently, 



m 



27ih 2 (3 



(57) 



with V(x) and U(q) defined in ( f40"|) and (|4l). This is, clearly, the "classical" limit for the 
partition function with a pre-factor that incorporates quantum fluctuations. 

At low temperatures, — > oo and ( |55| ) is solved for fee — > 1. A careful derivation |I[ 
leads to 



Z 2 (0) e ~°° 



9© rfg, 



' e -/bc]/s 



(58) 



with g@ given by 



1k'\ 



\l-2k' 



4V2e" / 2 , 



(59) 



and J[g c ] given by 



1 2 2 \ 



3/2 



+ O(0e- e ). 



(60) 
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3. Applications 



We shall now apply the quadratic semiclassical approximation to obtain the ground- 
state energy and the curve for the specific heat as a function of temperature. These two 
applications will teach us about the power of the approximation. 

In order to compare fl58|) with the expected low-temperature limit of the partition func- 
tion, Z(Q) ~ e _eeo ^ (where e (g) = E (g)/huj is the dimensionless ground state energy), 
it is convenient to rewrite it in a form in which the G-dependence can be analyzed more 
easily. This can be done by changing the integration variable back to q . Since q t nctie = qo 
and g@ is the value of q t corresponding to q — > oo, one has 



Z 2 (8) 6 ~°° 





exp < — — 



1 \ 3 / 2 
1 + 7^ ' 



(61) 



When ^> 1 it is possible to write an approximate expression for qt(qo, @) JTJ] , thus allowing 
to write the integrand in (|B~TD solely in terms of go an d O. The final result is 



Z 2 (G) 



e^oo 



2e"©/2 roo 
JlU) J-oc 



dq 



exp < 


r 4 






3/2 

-1 


} 


v /l + l% 2 


(l + V 1 + 5 9?) 





(62) 



This gives Eo(g) = 1/2, indicating that the quadratic approximation is insufficient to yield 
corrections to the ground state energy of the harmonic oscillator. On the other hand, if one 
recalls that the partition function can be written as 



z(e) 



p(Q;q,q) dq, 



(63) 



where 



p(©;g,g) = E e ~ ee " \Mq)\ 



2 e 



oo 

~ e 



l^b(?)| s 



(64) 



is the diagonal element of the density matrix, one may take the square root of the integrand 
in (|62D as an approximation to the (unnormalized) wave function of the ground state. To 
test the accuracy of this approximation, we have evaluated the expectation values of the 
energy for some values of g and compared them with high precision results found in the 
literature. As Table | shows, the ground state energy computed with this "semiclassical" 
wave function differs from the exact one by less than 1% even for g as large as 2. 

Another concrete problem that can be treated is the calculation of the specific heat of 
the quantum anharmonic oscillator. It can be written in terms of Z(Q) as 



i d 2 z _n dzx 

~Zd& I ~\ZdQ) 



(65) 



This expression was computed using MAPLE for a few values of G and the coupling constant 
value g = 0.3. The result is depicted in Fig. H, which also exhibits the curve of specific heat 
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of the classical anharmonic oscillator (solid line). As expected, the results agree when the 
temperature is sufficiently high, but, in contrast to the classical result, the semiclassical 
approximation is qualitatively correct at low temperatures too, dropping to zero as T — > 0. 

This result, together with the estimate for the ground-state obtained previously, shows 
that the quadratic approximation works very well, being quite accurate at high temperatures, 
and still reliable at lower temperatures. In the next subsection, we will comment on why 
this is so. 



4- Beyond quadratic 

In this subsection, we shall compute a first correction G\ to the quadratic approximation, 
which corresponds to the m = 1 term in flT2"p. Using Q5DD, we obtain 



Z(Q) = Z 2 (S) - - / dq e- I[qc]/9 / dO 
g J-oo Jo 



+ .... (66) 



Eq. (H) yields (r? 3 ) = and (r/ 4 ) = 3g 2 G c (0, 0; 9, 0) G 2 (9, 9). Inserting these results in (|| 
and changing the integration variable from go to qt gives 



Z(&) = r dq t D(q u 9) e- 1 ^ [1 - ga x {q t , 9) + . . .] , (67) 

where 

ai(g t ,9) = ^ [ d9g 2 (9,6). (68) 



A Jo 

Because of the complicated form of Q(9,9), it is not a simple task to compute ax(q t , 9). 
However, we can estimate the magnitude of this term without much effort. Indeed, as 
shown in |l|], Q(9,9) obeys the following inequality: 

G(9,9) < 6{e ~ e) (0<9<Q). (69) 

Therefore, 

oi( ft) e)<^. (ro) 

This shows that this correction to the quadratic approximation, Eq. (|51|) or (|53|), can be 
neglected whenever the condition g9 3 /40 <C 1 is satisfied; this is compatible with the 
numerical agreements obtained in the applications of the quadratic approximation. 

The next term in the expansion for Z(Q), which corresponds to the m = 2 term in (|12|). 
has a piece with a factor g and one with a factor g 2 . The former comes from the product 
of (t] 6 ) ~ g 3 with the overall g~ 2 , whereas the latter involves (^ 8 ) ~ g A . The Feynman 
diagrams which correspond to the m = 1 and m = 2 contributions are depicted in Fig. |j. 
Note that two of the m = 2 diagrams involves the three- leg vertex in floTf), which depends 
explicitly on q c {9). 9rdinary perturbation theory corresponds to the subset of graphs which 
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do not contain the three-leg vertex with the replacement of Q (9, 9') by the corresponding 
(free) expression for the harmonic oscillator. In fact, Q{9, 9') can be expanded in terms of its 
(free) harmonic oscillator expression and insertions of (U"[q c ] — 1), already an indication of 
its non-perturbative nature. Alternatively, we may obtain the perturbation theory diagrams 
by letting q c {9) — > 0. 



IV. CONCLUSIONS 



The results of section [TJ can be generalized to higher- dimensional Quantum Statistical 
Mechanics, just as in Quantum Mechanics, where this was accomplished in ||T0| , [TT|j . The 
generalization to potentials which allow for more than one classical solution, such as the 
double-well quartic anharmonic oscillator, requires a subtle matching of the series around 



each appropriate saddle-point (i.e., the minima). This is presently under investigation [£1 



An extension of our results to field theories is hampered by the fact that we do not know 
how to construct semiclassical propagators in general. The technical simplifications which 
appear in Quantum Mechanics cease to exist. However, our methods may still be of use in 
problems where classical solutions have a lot of symmetry (e.g., spherical symmetry) so that 
we can reduce them to effective one- dimensional problems. There are many such examples 
in Physics: instantons, monopoles, vortices and solitons are a few of the backgrounds that 
fall into that category. We are currently pursuing this line of investigation. 

Finally, we should remark that the field-theoretic treatment can be used to compute 
any correlation function of interest, in the usual manner. Therefore, a semiclassical series 
can be written down for any physical quantity once it is expressed in terms of the relevant 
correlation functions. 
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APPENDIX A: 



It remains to show how one can obtain G c (9i, rji, 92, 172) from the classical path. For 
this, we use the fact that the action I2 is quadratic in 77, and so the path integral in (|20|) is 
completely determined by the extremum rj e (9) of /2[^i, #2! v]i which satisfies Eq. (p4|), subject 
to the boundary conditions r}(9i) = % and 77(^2) = Vz- Thus, 

G c (9 1 ,Vu02,m) = G c (9 1 ,0;9 2 ,0)e- I ^ e ^, (Al) 

where, after an integration by parts, 

h[9i, 9 2 ; Ve ] = 1 [ m fjM - Vi UOi)} ■ (A2) 
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We can obtain rj e (9) by finding the linear combination of any two linearly independent 
solutions, Tj a {0) and r]b(9), of (|24]) which satisfies rj e (9i) = 771 and rj e (9 2 ) = r] 2 . The result is 



with Q(9, 9') as defined in We may then write 

1 



-M^iA^e 



2fi 
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[Wia ^ + W21 77^ - (Wn + W 22 ) 77x772], 



(A3) 



(A4) 



where fiy = Q(9i,9j) and = dQij/d9j. (Note that Wa is the Wronskian of i] a and 77;, 
computed at 

Explicit expressions for ^(6*) and 776(6*) can be obtained as follows. By differentiating @ 
with respect to 9, one can verify that rj a {9) = q c {9) satisfies fl24"|). For the second solution, 
we take rj b (9) = q c {9) Q(9), where Q{9) is defined as 

r e d9' 

Q(9) = Q(0) + 



q. 



(A5) 



for 9 < 0/2, Q{9) = -Q(0 - 9) for 9 > 0/2, and Q(0) is chosen so as to make rj b (6) 
continuous at 9 = 0/2. One can easily check, using (^), that rjb(9) indeed satisfies (p4|). 
(Alternatively, one could use a procedure introduced by Cauchy [fT0| , |lTj| , and differentiate 
the classical solution q c (9) with respect to any two parameters related to its two constants 
of integration.) We can now write explicit expressions for Q12 and Wif 

ftia = UOi) U9 2 ) \Q{9 2 ) - Q(9 1 )}, (A6) 



Wi, 



qmU'[q c {9 3 )] [Q{9 3 )-Q{9 1 )} + 



q c 



(A7) 



c\ u 3> 



As a final step, the pre-factor in (|A1|) can be derived JTJ using the methods of Refs. 
23|,|20|. The result is 

1/2 



Gc(0i,O;0 2 ,O) 
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2?r#fi 
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(A8) 



From (|A7|), one easily finds = 1. Therefore, our quadratic semiclassical propagator is 
given by 



g c (0i, #2,772) 



exp 



2#ft 



(Wia 772 + ^21^- 2 77x772) 
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As promised, it is completely determined by the classical solution. 

Finally, we note that the van Vleck determinant A is a by-product of 

A(q Q , 0) = G c - 2 (0, 0; 0, 0) = 2ng £1(0, 0) = 4ng q 2 c (0) Q(0). 

As shown explicitly in one can express A as 

a ^g[U{q t )-U{q )\ fdG\ 



U'(q t 



dq t 



(A9) 



(A10) 



(All) 



Together with (||), this shows that one does not need to know q c {9) in order to write the 
first term in the semiclassical series; it is enough to know qt(qo, ©)• 
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TABLES 



TABLE I. Ground state energies for different values of g (h = m = u> = 1). 



g E^semiclassical) 3, E (exact) h error(%) 

04 0.559258 0.559146 0.02 

1.2 0.639765 0.637992 0.28 

2.0 0.701429 0.696176 0.75 

4.0 0.823078 0.803771 2.40 

8.0 1.011928 0.951568 6.34 



a (4>o\H\4>o)/((j)o\(j)o}, where (j)o(qo) is the square root of the integrand in Eq. fl62]). 
b Values quoted from Ref. [pSjj . 
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FIGURES 
FIG. 1. U{q). 

FIG. 2. Graph of f(k). 

FIG. 3. Specific heat vs. temperature (T = 1/0) for the quantum (diamonds) and classical 
(solid line) anharmonic oscillator, g = 0.3. 

FIG. 4. Feynman graphs for m = 1 and m = 2. 
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